Destructive quantum interference in spin tunneling problems 
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In some spin tunneling problems, there are several different but symmetry-related 
tunneling paths that connect the same initial and final configurations. The topolog- 
ical phase factors of the corresponding tunneling amplitudes can lead to destructive 
interference between the different paths, so that the total tunneling amplitude is 
zero. In the study of tunneling between different ground state configurations of the 
Kagome-lattice quantum Heisenberg antiferromagnet, this occurs when the spin s is 
half-odd-integer. 
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The problem of calculating the rate at which a quantum spin system tunnels between 
its different low-energy states has been of interest in various different contexts The 
tunneling amplitude is usually calculated by setting up a coherent-spin-state path integral 
and analytically continuing to imaginary time (t —>■ —ir); the leading contribution can 
be found using the method of steepest descent 0. The phase of the tunneling amplitude 
depends on a topological phase 

In this Letter we point out that when there are several different symmetry-related tunnel- 
ing paths connecting the same degenerate classical ground states, the topological phase can 
lead to destructive quantum interference between their amplitudes, and hence to a total tun- 
neling amplitude that is exactly zero. The occurrence of such a cancellation has a geometric 
interpretation and depends on the value of the spin s. We also present two examples: The 
first involves a Hamiltonian with an m-fold symmetry axis. The second concerns tunneling 
amplitudes between different degenerate ground state configurations of a Kagome-lattice 
quantum Heisenberg antiferromagnet; these amplitudes are zero if s is half-odd-integer, but 
non-zero if s is integer. 

For extended quantum spin systems in one and two dimensions, qualitative differences 
between integer and half-odd-integer spin s have been reported previously [Q. In these 
papers, quantum interference between topologically distinct paths of a continuous unit vector 
field n{x,t) leads to interesting s-dependencies. Our examples show that similar effects occur 
for much simpler systems that involve only a small number of individual spins. 

The tunneling amplitude: In tunneling problems, the customary object of study is the 
imaginary time transition amplitude from an initial state \i) to a final state |/). For a spin 
system, this can be written as a coherent-spin-state path integral 0: 

Uf, = (/|e-^*/^K) = J Vne-^/"" . (0.1) 

where S = J drCis the Euclidean action andVQ is the measure of the path integral. For the 
special case of a single spin, the Euclidean Lagrangian is 

C = -ihs^{l - cos 9) + ni^, 9) . (0.2) 

2 



The coordinates {4>, 9) label the coherent spin state |0, 9) for a particle with spin s, and may 
be associated with a unit vector n in the (0, 9) = 0. direction. The dot on means dr- The 
"semi-classical" Hamiltonian is the expectation value 7i = {(l),9\'H\4>,9) of the operator H. 
For simplicity, we shall consider only the case where \i) and |/) are "classically degenerate 
ground states" (in the sense that (i|7i|i) = {f\H\f) are the smallest possible expectation 
values oiH), and are separated by an energy barrier. 

The above path integral can be evaluated by the method of steepest descent: 

Uf^^Y.-^^^-'^'^'^Euih (0.3) 

where H^^^ fD{n- n«)e-(^^^"^+^^^'"-)/^ . (0.4) 

Here S'^'> is the action evaluated along the l-ih "tunnehng path" , which is a solution to the La- 
grangian equations of motion and will be denoted by overlined variables, e.g. (0^'^ (r) , 9^^^ (r) ) . 
The index / allows for the possibility of several different symmetry-related tunneling paths. 
The prefactors AA*^'-* measure the effects of fluctuations around the Z-th tunneling path. 
Tunneling problems are characterized by the fact that the coordinates in general acquire 
imaginary parts along the tunneling path (else it is not possible to satisfy the requirement 
that the Hamiltonian be a conserved quantity along the path). Consequently, the vari- 
ous S'p can have non-zero imaginary parts. Quantum interference, and possibly complete 
cancellation (so that Ufi = 0), can thus occur between the amplitudes of the different paths. 

Geometric interpretation of phase: The conserved energy Eg — 7^(0*^^'*) along the tun- 
neling paths that connect the degenerate states \i) and |/) may be set equal to zero without 
loss of generality. Hence 5^'-' is completely determined by the first term of C This term has 
a well-known geometrical interpretation, which we now discuss. 

Let ^{t) — {(I){t),9{t)) be a purely real, closed, non-self-intersecting, smooth path in 
spin space. The area on the unit sphere enclosed by this path is given by 

A = JdT^{l-cos9) , (0.5) 
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modulo 4:71, depending on which of the two oriented areas (the "inside" or the "outside" of 
the closed path) one considers. 

Now consider two (of the possibly many) tunneling paths, say fi^''*(r) and il^^ ^(r), not 
necessarily purely real, connecting \i) and |/). Suppose that some symmetry of the Hamil- 
tonian ensures that the absolute values of the two tunneling amplitudes uf- and t/j-'j ■* of 



eq. (pTSl) are equal, i.e. that 

Re[Si'^] = Re[SP] , and AT^') = Af^''^ . (0.6) 

Intuitively speaking, this situation arises when the local neighborhoods of the two tunneling 
paths are, for symmetry reasons, identical for the two paths, so that the local shapes and 
sizes of the barriers (which determine the Re[5^'^]'s) and the local fluctuations around the 
tunneling paths (which determine the A^*-''''s), are identical. Then f/j'^ and U^''^ differ at 
most by a phase, namely 



{Si'^ - SP)/h = -IS JdrRel^^'^] (l - Reicos^^'^] 
fdrRelP'^] (l-Re[cos^('') 



(0.7) 

Terms such as Re[ ]Im[ ] do not appear due to eq. ( p.6| ). We have assumed that 
Im[0(')]Im[cos5W] = (same for / <-> Z'), for the following reason: The conserved energy 
condition, H(n^'V) ) = 0, can be solved to find, for example, 6^''^ in terms of (f)^''\ The 
dependent variable O^''^ will be a complex function of the independent variable (f)^^\ which 
can be taken to be real, so that lm[0*^'''] = 0. A similar argument works if one chooses to 
write 0*'''' in terms of 6'*'''' [|]. 

Now, let All' be either one of the two oriented areas on the unit sphere bounded by the 
loop formed by the two paths Re[f2'"'''] and Re[f2'^' ■*]. Then, by eq. (p.5|) , the relative phase 
eq. (pT^) reduces to [H 



{Si'^-SP)/h = -tsAi, . (0.8) 
The 47r-ambiguity in Aui is irrelevant, since exp(— z47rs) = 1 for any spin s. 



If sAiii is an odd multiple of tt, the amplitudes for the two paths interfere destructively, 
so that Uj^ + Uj\ ^ = 0. In the simplest case where fi*'''' and fi^' ^ are the only two tunneling 
paths, this means that the total amplitude Ufi is zero. Note that this result does not depend 
on the detailed dynamics of the tunneling motion. In particular, it is not necessary to know 
the imaginary parts of the tunneling paths fi''''' explicitly, as long as symmetry arguments 
can be invoked to assert that the relations ( p.6|) hold. It may happen that symmetry ensures 
that the prefactors A/" are equal to all orders of the steepest descent method, which is an 
expansion in powers of 1/s (and not merely to the lowest order that is usually employed). 
In this case the cancellation, if it happens, is exact to all orders in 

Clearly, the occurrence of destructive interference depends crucially on the value of s. 
This could create complications in the analysis of the tunneling behavior of, for example, a 
small ferromagnetic particle with effective spin ns (as has been considered in c]), where 
n is the number of spins in the particle, since the effective spin is then n-dependent. 

Turning on an external magnetic field can in general affect the occurrence of destructive 
interference by changing the initial conditions (i.e. \i) and |/)) and the area Aai enclosed 
between different tunneling paths, and by destroying the symmetry that ensures that the 
absolute values of the amplitudes for all tunneling paths are the same. 

We now present two examples where destructive interference does occur. 

Hamiltonian with m-fold symmetry axis: Consider a single spin s. Suppose that the 
Hamiltonian has an easy axis (say the ;z-axis, 6 = 0), around which it has m-fold rotational 
symmetry: H{(p, 6) = 7-^(0 + 27r/m, 6) for all 0. Suppose that \i) = \9 = 0) and |/) = \9 = n) 
correspond to the two degenerate classical ground state configurations. Clearly, if fi*"^'' = 
{(j)^^\6^^^) is a tunneling path from \i) to |/), so are the paths fi''''' = (0*'°'' + 2nl/m,9^^^), 
for / = 0, . . . , m — 1 (see fig. |T|). By symmetry, all Re[5y^] are equal to each other, as are all 
A/''-'''. Furthermore, 

(5(°) - Sll^)/n = -isAoi = -isAnl/m , (0.9) 

because the area on the unit sphere enclosed between the real parts of any two neighboring 
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paths is, by symmetry, necessarily equal to An/m. The total amplitude is thus: 



m— 1 



1=0 



where the 5-function is non-zero only if 2s is an integer multiple of m. 

A realization of the above scenario, with n = 3, is afforded by a Hamiltonian with cubic 
anisotropy and a small additional uniaxial anisotropy along the [111] direction (thus 3- fold 
symmetry around [111]): 



where n^, Uy and are the direction cosines of the unit vector n in the (</>, 9) direction, and 
Ji ^ J2 > 0. The unit vectors d = :^(1, 1, 1) and —d define the two classical ground-state 
configurations 

Spin tunneling in the Kagome lattice: Our second example concerns spin tunneling events 



where srij = 6'j|Sj|0j, 6'j) is the "classical" spin (see fig. |^a). Any configuration in which 
the spins on each triangle minimize their energy by assuming a coplanar configuration, with 
relative angles of 120°, (see fig. §b), is a classical ground state. Therefore there are macro- 
scopically many degenerate classical ground states. Generally fluctuations, both quantum 
and thermal, are expected to lift some of the degeneracies, thereby inducing magnetic order- 
ing by selecting particular configurations ("order from disorder"). For example, spin wave 
expansions about various ground state configurations have shown that maximally copla- 
nar configurations, in which all spins in the lattice lie in the same plane (let this "reference 
plane" define = and (p = n), have lower zero-point energies than any other configurations 



On the other hand, tunneling between different degenerate ground state configurations 
competes with "order from disorder" selection effects, because it tends to drive the system 




(0.10) 



in a 2D quantum Heisenberg nearest-neighbor antiferromagnet on a Kagome lattice (fig. 0) 
1^. The Hamiltonian is taken to be 

7i = sV^ni-nj, (J>0), (0.11) 



0, b,d]. 
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into a superposition of degenerate states, rather than selecting a particular one. As the 
simplest example of a tunneling event on the Kagome lattice, we consider the so-called 
"weathervane mode" (see fig. |^): the six spins of an AB AB AB-hexa.gon in one maximally 
coplanar ground state, rotate synchronously by 180° around the z-axis (defined by the 
C-spins), to end up as a BABABA-hexagon in another maximally coplanar ground state, 
I/), while all other spins remain fixed. Due to the above-mentioned spin wave selection 
effects, there is a "coplanarity barrier", say Jb{<f)), to this type of motion, which, by reflection 
symmetry in the reference plane, obeys Jb{<p) = Jb{—(p)- 

To study the hexagon tunneling event, we consider the following Euclidean Lagrangian: 

6 _ 6 

C = '^—ihs(j)j{l — cosOj) + sin^i(0j + (f)j^i + tt) 

6 

+ s^J^[nj-nj+i + nj-z + rij+i-z + |]. (0.12) 

The index j is defined modulo 6. We take J ^ Jb > 0. The s^Jrij-z terms (interaction 
with C-spins) and the phenomenological sin^0 coplanarity barrier are assumed to be the 
only ways the other spins in the lattice, which are assumed to remain fixed, influence the 
six spins on the hexagon. 

To minimize the cost of the dominant s^J term, the hexagon spins are expected to 
rotate collectively, maintaining mutual near-coplanarity. Indeed, it can be shown [0 that 
the hexagon tunneling problem can be mapped onto a simple model problem, defined by the 
Lagrangian 

£= -ms4){l - cos 9) + 12sV(cos 9 - \ f + 6Jb sin^ , 

involving only a single (collective) spin degree of freedom with an effective spin of 6s. For 
present purposes, however, the following observations suffice (discussed in detail in [0): 
Due to the reflection symmetry about the reference plane, there are two possible tunneling 
paths, to be denoted by (0^,F^); they differ from each other only in the direction of the 
0-rotations and satisfy 0^- (r) = — 0^(r), for j = 1, . . . , 6. Reflection symmetry ensures that 
Re[5^] = Re[5;r] and J\f^ = J\f^ . Along both tunneling paths, every (p^ is purely real. 
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To satisfy H = during the tunnehng event, each 6j develops a time-dependent imaginary 
part (which vanishes in \i) and |/)), but Re[cos6'j] maintains the value it has in \i) and |/), 
namely Re[cos 9j] = — ^ j = 1, . . . , 6. 

Thus, for each of the six spins, the real part of the tunneling path is a contour of constant 
Re[cos6'j] = — |, with Re[0j'^] (or Re[0j ]) changing from to tt (or — vr) for (+)- or (— )- 
paths. For each spin, the area enclosed between the (+)- and (— )-paths is thus equal to vr, 
giving for the six spins a total phase difference of iGns between the amplitudes for an (+)- 
or (— )-event. It follows that the total tunneling amplitude becomes 

Ufi = Ar+e-'5°^/^(l + e'^"'^) (0.13) 
2J\f ^ e~^^ if s is integer , 
if s is half-odd-integer . 

Similarly, consider any larger closed "loop" of alternating A- and i?-spins within a 
ground-state configuration. It can be proven that any such loop contains 4n + 2 spins 
(n is some integer) |jTO|. Again one can study the tunneling between two configurations that 
only differ by 0; — ^ 0; -|- vr (i.e. A ^ B) for each spin on the loop. The relevant phase 
between (+)- and (— )- paths will be m{An + 2)s, and for half-odd-integer s, destructive 
interference again occurs. 

The above results imply that in that subset of parameter space where "order from dis- 
order" selection effects and the competing tunneling effects that favor more disorder are 
more or less equally important, one might expect interesting integer versus half-odd-integer 
s effects, reminiscent of those found in ID antiferromagnetic spin chains [^, a]. It would 
therefore be of interest to find integer-spin realizations of the Kagome lattice antiferromag- 
net, since s is half-odd-integer for the two experimental realizations of the Kagome lattice 
that have been proposed: the magnetoplumbite SrCr^^xGca+xOi^ (s = |) and the second 



layer of He^ atoms on a graphite substrate (s = |) [|TI|. 

Integer versus half-odd-integer s effects might also make their appearance in exact diag- 
onalization studies of finite size systems with discrete degeneracies, in the analysis of which 
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spin tunneling methods should be very useful. The simplest realization of such a discrete 
degeneracy is the J1 — J2 square lattice antiferromagnet for J2 > Ji/^, which shows a discrete 
degeneracy between antiferromagnetic ordering vectors (1,0) and (0, 1) ||12|| . 

The vanishing of tunneling amplitudes obviously implies an exact ground state degen- 
eracy in the semi-classical picture. Sometimes, this degeneracy can be shown to exist for 
all eigenstates of the system, on a purely quantum-mechanical level. For example, consider 
the following toy model for a quantum spin (independently suggested to us by V. Elser): 
Ti. = —Sl — aS"^. Here a -C 1, and Sz and Sx are spin operators. The tunneling am- 
plitude between the two classical ground states, n = ±z along the two tunneling paths 
6'^ : ^ TT, 0''' = 0, = TT, is zero for half-odd-integer s (the phase difference is 27rs). 
On the other hand, Ti has time-reversal symmetry and hence all its eigenstates display a 
two-fold Kramers degeneracy for half-odd-integer s. It would be interesting to investigate 
more generally under what circumstances vanishing semi-classical tunneling amplitudes also 
imply exact quantum-mechanical degeneracies for all eigenstates. 

In conclusion, we have shown that the topological phase factor occurring in spin tun- 
neling amplitudes can have quite striking effects. If there are several symmetry-related 
tunneling paths whose individual tunneling amplitudes have the same absolute value, their 
topological phases can lead to destructive quantum interference between the paths and a 
total tunneling amplitude that is zero. The conditions under which this occurs can be in- 
terpreted geometrically in terms of the areas on the unit sphere enclosed between the real 
parts of the various tunneling paths. In quantum antiferromagnets, interesting integer ver- 
sus half-odd-integer s effects can result from the competition between "order from disorder" 
selection and tunneling effects. 

Helpful discussions with A. Auerbach, D.P. DiVincenzo, V. Elser, A. Garg, M. Kvale 
and S. Sachdev are gratefully acknowledged. Related preliminary calculations were done by 
Q. Sheng. Part of this work was supported by N.S.F. grant DMR-9045787. 
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FIGURES 

FIG. 1. Three equivalent tunneling paths for a Hamiltonian with 3-fold rotation symmetry on 
a "wrapped-open" unit sphere. 



FIG. 2. An ABABAB hexagon of spins (on the left) in a coplanar ground state configuration. 



FIG. 3. (a) The unit vector n = {(f), 6). (b) Type A, B and C spins on a triangle of the Kagome 
lattice. 
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